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The Deutsch-Jozsa algorithm is experimentally demonstrated for three-qubit functions using pure 

coherent superpositions of Li2 rovibrational eigenstates. The function's character, either constant 

or balanced, is evaluated by first imprinting the function, using a phase-shaped femtosecond pulse, 

on a coherent superposition of the molecular states, and then projecting the superposition onto an 

- - - ionic final state, using a second femtosecond pulse at a specific time delay. 
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Quantum computation R-fl] is aimed at utilizing the quantum nature of physical systems in order to solve com- 
putational problems in efficient ways that are impossible in classical computation. One of the benchmark quantum 
algorithms is the Deutsch-Jozsa (DJ) algorithm M. Its task is to distinguish whether a binary n-qubit function 
Ps) / : {0, 1}" -^ {0, 1} {n is given) is constant or balanced. The output of a constant function is identical for all possible 

^ n-qubit inputs, while the output of a balanced function is 1 for half the possible inputs and for the other half. 
OO The DJ quantum algorithm identifies the constant/balanced character of a function (known to be either constant or 
balanced) in a single call to the function, as compared to the corresponding classical algorithm that requires 2"^^+l 
function evaluations to provide a solution. This improvement results from the inherent parallelism when applying 



o 

^— V ' a function-dependent unitary transformation on all the possible input elements, which are simultaneously contained 
^_^ within a coherent superposition. 

f^ In recent years, the Deutsch-Jozsa algorithm, in its revisited form B, has been implemented experimentally mainly 

^5 , with NMR techniques. Using pseudopure states |q], it was demonstrated for functions having up to a four qubit 

Q," input 0]. Also, an implementation of the algorithm for two-qubit functions using single photon linear optics has 

I ' been published [p[. A modified version of the algorithm, which does not require a control qubit, has been proposed 

rH p and implemented using NMR for functions with inputs composed of up to three qubits Q- In this version of 

d the algorithm the coherent superposition representing the function domain is obtained by applying a first Hadamard 

^ rotation on the initial state |0)|0) • • • |0) of the n-qubit system. The unitary transformation constructed for each of the 

^ '' various functions introduces a function-dependent phase to each of the elements of the superposition. After a second 

^ , n-qubit Hadamard transform is applied, the superposition either ends in the initial state if the function is constant 

or in any other state if it is balanced. In the present work, this modified algorithm is implemented for three-qubit 

functions using pure coherent superpositions of molecular eigenstates (wave packets) . 

jrt Our main motivation behind the current work is exploring the use of a new experimental platform for small-scale 

" " ' problem-specific implementation of quantum algorithms. The method is based on an ensemble of small molecules in the 

gas phase interacting with a sequence of multiple shaped femtosecond laser pulses. The computational task is carried 

out through the time-dependent dynamics of the molecule, with short computation time and very low decoherence rate. 

The system is characterized by a value of lO'^-lO'* for the ratio between the decoherence time and the manipulation 

time by a shaped femtosecond laser pulse. The experimental technique is such that the computation involves pure 

states only, a characteristic that is highly important for a physical implementation of quantum computation. These 

pure molecular states span several (entangled) internal degrees of freedom - rotational, vibrational, and electronic. 
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Currently, we aim at providing a physical platform for a problem- specific quantum computation. This is different 
from the common formalism of a universal quantum computation [^ , which is also the one most physical implementa- 
tions of quantum computation have followed. Under the universal formalism, every unitary transformation acting on n 
qubits is decomposed into a quantum circuit, which is a sequence of elementary gates, each acting on one or two qubits. 
The gates composing the circuit are all members within a universal finite set of basic gates. This model provides 
a powerful universal programming language for quantum computation. However, it acquires a significant drawback, 
since most unitary transformations can only be implemented inefhciently, i.e., they require a circuit of elementary 
universal gates whose size is exponential in n. Only special transformations can be decomposed into polynomial-size 
circuits. One prominent transformation is the quantum Fourier transform, which is also a basis for the most successful 
quantum algorithm known today, that is the Shor's quantum algorithm for the factorization of numbers |ll| . Thus, 
it seems highly desirable that, in parallel to efforts along the above universal quantum computation model, other 
frameworks for quantum computation will be explored, even if they are non-universal. Our objective is to specifically 
design and shape laser pulses such that under the action of a single shaped pulse the molecules will undergo a complete 
specific unitary transformation. This might allow, at a price of losing the universality, the implementation of specific 
quantum algorithms in a much more efficient way as compared to using other physical platforms designed for universal 
quantum computation. The potential for achieving a complete transformation using a single pulse (or a very small 
number of pulses) lies in the richness and complexity of the interaction of molecules with a strong broad-bandwidth 
laser field. The corresponding dynamics involves simultaneously single- and multi-photon, direct and Raman, res- 
onant and off-resonant transitions among a large manifold of quantum molecular states. It is expected that useful 
theoretical computational methods for designing such pulses will include also optimal control theory |12]. Initial 
studies in this direction point to its feasibility [|3|jlJ] . The main limitation of the new platform as implemented in 
this work is the limited size of the experimentally accessible part of the whole Hilbert space of the internal molecular 
degrees of freedom. Current experimental techniques provide a basis for handling computational tasks involving pure 
states with up to hundreds of eigenstates (the equivalent of a Hilbert space of 8-10 qubits) spanning several internal 
degrees of freedom. This by itself is significant for the experimental research of quantum computation. Moreover, 
a potential for larger scale quantum computation lies in developing more sophisticated molecular excitation schemes 
and /or combining the present method with other techniques involving, for example, the trapping of cold molecules 



L5|,|l6| or cold molecular ions |17|. 



II. ALGORITHM IMPLEMENTATION 
A. Using Coherent Molecular Superpositions 

In the current implementation of the Dcutsch-Jozsa algorithm, 2" rovibrational eigenstates are being used to 
represent the eigenstates of a system of n entangled qubits. Each n-qubit function is represented by a different unitary 
operation stored in an oracle, and the task is to determine whether the unitary operation applied to the molecular 
system corresponds to a constant or balanced function. The various unitary operations are implemented using 
phase-shaped broadband femtosecond pulses. The interaction of the shaped pulse with a molecule, which is initially 
prepared in a pure single rovibrational state, transforms the molecule into a corresponding coherent superposition of 
rovibrational states. Different functions result in different superpositions. The physical features of the molecule-field 
interaction allow us to apply the first Hadamard rotation and the subsequent unitary transform, representing one of 
the functions, in a single step. The former is associated with the amplitude transfer and the latter corresponds to 
the phase modification: UphsUamp\tpit = 0)}. The evaluation of the function's character is accomplished by probing 
the overall molecular superposition through its projection onto an ionic final state. This is accomplished by ionizing 
the molecule at a single given time delay after excitation using a second (unmodified) femtosecond pulse. The ability 
to probe the overall superposition directly as a whole at one time, which originates from the quantum nature of the 
process, allows the readout step of the algorithm to be achieved without applying the second Hadamard rotation. 
The algorithm is performed on an ensemble of molecules. Overall, only one (encoding/imprinting) unitary operation 
and one measurement suffice to find out the function's character. 

Although, in general, 2" molecular levels are sufficient to represent the Hilbert space of n qubits, in the present 
experiment 2" -I- 1 states were used. They include 2" rovibrational levels in the electronic E-state of Li2 and a single 
initial rovibrational level in the electronic A-state (see Fig. |]). This simplifies the implementation of the unitary 
operations without any principal limitation. Here, the electronic correlation between the E-state and the A-state 
is not probed |lq], and, thus, the Hilbert space naturally reduces to the 2" levels in the E-state. The algorithm 
is demonstrated here with three-qubit functions, n = 3, thus eight rovibrational states are employed in the E- 
state. Each rovibrational level, (ve, Je), corresponds to an eigenstate of the threc-qubit product space, denoted as 



|fc) = |77i)|n)|o) with m,n,o = {0,1} where k is the decadic representation of a three-bit binary digit. ExpHcitly, 
this means that in the present experiment \ve = 13, Je = 17) = |0) = |000), \ve = 13, Js = 19) = |1) = |001), 
\VE = U,Je = 17) = |2) = |010), \vE = U,Je = 19) = |3) = |011), \ve = 15, Je = 17) = |4) = |100), 
\VE = 15, Je = 19) = |5) = |101), \vE = 16, Je = 17) = |6) = |110), and \ve = 16, Je = 19) = |7) = |111). 

The unitary transformation, which represents a constant/balanced function, is induced in the molecule by the 
shaped laser pulse via the Hamiltonian 
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The symbols fifc = fik£k{t)exp{~iLUkt) (k = — 7) are the amplitudes acquired by the rovibrational levels on the E-state 
following their excitation by a transform limited (i.e., phase- unshaped) pump pulse from the initial rovibrational level 
on the A-state. Each amplitude flk is determined by the transition dipole moment fi^ between the initial state and the 
excited state fc, and by the spectral magnitude efc(t) of the field at the specific excitation frequency ujk between these 
two states. The time-dependence of the spectral magnitude reflects the pulsed character of the field. The quantities 
Cfc = exp{—i(f)k) include additional phase factors, 0^, that are introduced to the excited rovibrational level by the 
phase shaping of the pump pulse. In the current description the influence of off-resonant coupling is neglected. The 
Hek are the field- free Hamiltonians of the various excited rovibrational levels, and they arc equal to the energies of 
the levels, i.e., Hek — hujk- Similarly, Hg corresponds to the initial level on the A-state. Under weak field conditions, 
the excited wave function on the electronic E potential can be formulated by first order perturbation theory [|l8|-p0[ . 
After eliminating the A-state dynamics by putting Ha — 0, the excited superposition on the E-state is given at time 
T after excitation, when the pump pulse is over, as |lq| 
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As noted above, k = {ve,Je) while \k) represents the rovibrational eigenfunction of state k on the E-state. The 
molecule-field interaction excites each molecular level with a phase and amplitude that are controlled experimentally 
by shaping the excitation pulse. 

Following the formulation of the modified Deutsch- Jozsa algorithm 0] , the present pulse shaping is carried out to 
correspond to a function / such that 
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where a^ denotes the complex amplitude for a level k on the electronic E-state, which is independent of the specific 
evaluated (constant/balanced) function. Note, not all the ak are necessarily equal. The term (— l)/(l'=)) introduces 
the function dependent phase factor. The and 1 values of the function are encoded as a phase of 0° or 180° (i.e., 
a +1 or —1 factor), respectively. As a result the expansion in first order perturbation theory, [/ = 1 — iHt, provides 
a model for arbitrary phase and amplitude transfer, U — Uphs Uamp, from the initial A-state level onto the E-state 
rovibrational wave packet. It can be viewed as follows: 
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The right term of the r.h.s., Uamp, represents the function independent part, while the left term on the r.h.s., Uphs, is 
the unitary operation that encodes the function values as phases. The use of perturbation theory does not represent 
any principle limitations, since the theory of coherent control, for a closed quantum system of a discrete spectrum, 
ensures that complete control can always be achieved even with strong fields ||2^ . 

Each three-qubit binary function is given as a set of eight binary (0 or 1) values, each corresponding to a possible 
state, |fc) = \m)\n)\o), of a three-qubit input. There are 72 three-qubit constant/balanced functions - 2 constant and 70 



balanced. As mentioned above, a function's value of 1 is represented by a phase value of 180° and a function's value of 
by a phase value of 0. Following the oracle operation for a given function, a set of eight phases is determined. Those 
phases are then encoded experimentally into the pump pulse that excites the molecular superposition. In practice, in 
the current experiment, we have chosen to encode these phase values into the phase-shaped pump pulse as an increment 
over a basic initial set of phases apphed to the excited rovibrational states, <i>^°^ = {</'„J Je^ ^s = 13—16, Je = 17, 19}. 
This <i?^°^ set of phases is actually part of the function independent a^ coefRcients introduced above. As a result of 
this procedure, the two constant functions, /i = {0,0,0,0,0,0,0,0} and /2 = {1,1,1,1,1,1,1,1}, are represented 
by molecular superpositions having <I>(°) and ^'^"^ + 180°, respectively, as their 0^ set of phases (see above). Since 
the measured pump-probe ionization signal is sensitive only to the relative phases between the various wave packet 
components, i.e., insensitive to a global phase of the wave packet, both constant functions correspond to the same 
pump-probe signal. The specific set of phases *(°) used here is $(°) = {(fifln = 298.1°, (/.^"^ig = 352.0°, </)^4^i7 = 

215.9°, (/-J^'^ig = 137.9°, (/.^g^^iy = 169.7° , ^Jg^^ig = 337.6°, (/.igliy = 192.1°, ^ie.ig = 0°}. These values were chosen such 
that the ionization of the corresponding Li2 wave packet at 5 ps delay time after its excitation will result in a global 
maximum of the measured coherent signal. This is based on previous detailed coherent studies we have conducted on 
Li2 |iyj20]. The 5 ps time was chosen arbitrarily. The balanced functions (/s to 772) will result in rovibrational wave 
packets having sets of relative phases that are different from <i>(°\ such that the corresponding amplitudes of the wave 
packet ionization signals at 5 ps delay time will be significantly lower than the global maximum signal. Hence, the 
identification of the function's character can be made by measuring the signal amplitude at this single predetermined 
delay time, after a calibration control experiment to measure the global maximum signal (i.e., for the wave packet 
with $(*'') is initially performed one time only. 

B. Experimental Technique 

The molecular excitation scheme of the experiment, with the relevant potential energy curves of Li2 and Lij , is 
shown in Fig. H. The experiment |jl9,pQ| is conducted in a static cell, heated to 1023 K, that contains a lithium metal 
sample with Ar buffer gas at about 3.7 Torr (493.3 Pa). Using a cw single-frequency dye laser, an individual state-to- 
state transition, from a rovibrational level in the ground X^I]+ state to a rovibrational level in the excited A^S+ state, 
is induced on some of the thermally-populated ground state lithium dimers. As a result, a small fraction of the overall 
molecular ensemble is populating the single rovibrational state A^E+(wyi — 14, J a — 18) as a pure initial state. The 
information that represents a function is encoded into a phase-shaped femtosecond laser pulse using a pulse shaping 
setup. It employs a liquid crystal spatial light modulator (128 independent pixels) to control the spectral phase of 
the various frequency components of the pulse. The wavelength resolution of the pulse shaping setup is ~4.2 cm~^ 
per pixel, and the accuracy of a phase (0-360°) applied with the spatial hght modulator is better than 1°, depending 
on the particular phase value. The interaction of this shaped pulse with the Li2 molecules, populated in the selected 
initial [va = 14, J a = 18) state, results in a function-dependent tailored coherent superposition (wave packet) on the 
E^I]+ electronic state (pump step), which for three-qubit functions is composed of (ve = 13 — 16, Je = 17, 19) (eight 
states) . The wave packet excitation can be described within the weak field limit [Q . 

Using a second (unshaped) femtosecond pulse, this function-dependent excited wave packet is probed after a pre- 
selected delay time through the ionization of the molecule. The resulting pump-probe photoionization signal (ion and 
electron current) is the measured experimental quantity. It is measured using a pair of biased electrodes located inside 
the lithium cell and a lock-in detection system. All the laser beams (cw, pump pulse, and probe pulse) are periodically 
modulated using mechanical choppers. The measured signal contains a well-identified component that originates only 
from those molecules that have undergone excitation and ionization due to the timely ordered absorption of three 
photons: the first from the cw laser, the second from the pump laser pulse, and the third from the probe laser 
pulse. Overall, although most of the molecular ensemble is thermally populating the ground electronic state of Li2, 
the experimental technique allows an ionization signal to be obtained that originates only from those molecules that 
performed the quantum computing operation. The sub-ensemble containing these molecules is in a pure coherent 
quantum state. In general, the measured ionization signal is composed of a constant signal level and a part that 
depends on the pump-probe delay time |19| , pO[ . Currently, the measurement was averaged over 10^ pump-probe 
sequences at a specific pump-probe delay time. The pump and probe pulses originate from a Ti:sapphire laser system 
(200 kHz repetition rate) with ~160 fs duration, 12430 cm~^ central spectral frequency, ^^150 cm~^ bandwidth, 
parallel polarizations, and energies of ~0.5 and ~1.0 /xJ per pulse, respectively. 

The experimental conditions are such that the decoherence of the wave packet occurs on a time scale longer than 
5 ns. It is primarily due to collisions between the Li2 molecules and Ar and Li atoms, which result in pure dephasing 
and/or state-changing transitions. The decoherence time scale is at least three orders of magnitude longer than the 
excitation process that encodes the quantum information. 



III. EXPERIMENTAL RESULTS 



Figure displays several pump-probe ionization transients out of the complete set of measurements. They 
originate from various phase-tailored molecular wave packets, each corresponding to a different three-qubit con- 
stant/balanced function. For clarity, the transients in the figure are translated vertically to separate one tran- 
sient from the other. The baselines of all the measured transients, indicated in the figure by t he thin horizon- 



tal solid lines, are of the same value. Similar to the pair of constant functions (see Sec. II A ), following from 
the insensitivity of the measured ionization signal to a global phase of the wave packet, each transient actually 
corresponds to either of two functions, fi and fj, related as fi = fj, i.e., and 1 exchange in their logical rep- 
resentation (for example, /s and f^ given below). Such a relationship means that if fi corresponds to a set of 
phases $*, fj corresponds to ^^ = <I>* -I- 180°, i.e., an additional 180° global phase. The transient shown in all 
panels by the dashed lines corresponds to the constant functions (fi or {2), while the transients shown in thick 
solid lines each correspond to a different pair of balanced functions. The eight balanced functions presented in 
the figure are /a = {0, 0, 0, 0, 1, 1, 1, 1}, f^ = {1, 1, 1, 1, 0, 0, 0, 0}, /s = {0, 1, 1, 1, 1, 0, 0, 0}, /e = {1, 0, 0, 0, 0, 1, 1, 1}, 
fr = {1, 1, 0, 0, 0, 1, 1, 0}, /a = {0, 0, 1, 1, 1, 0, 0, 1}, /g = {1, 1, 1, 0, 1, 0, 0, 0}, and /lo = {0, 0, 0, 1, 0, 1, 1, 1}. The com- 
plete set of measurements contains results for all the 70 three-qubit balanced functions, measured at 5 ps delay time. 
In all cases, as can also be seen in the figure, measuring the coherent signal originating from the ionization of the 
wave packet at the single delay time of 5 ps provides the answer to whether the evaluated function is constant or 
balanced. The transients in the figure are shown over the full 2-8 ps delay time range only for clarity. The signal 
amplitude at 5 ps that corresponds to functions /a and f^ is the closest one, among all the 70 balanced functions, 
to the signal amplitude at 5 ps of the constant functions (/i and /2). Thus, overall, considering the experimental 
signal-to-noise levels with the corresponding error bars, the present probability of obtaining a correct answer about the 
constant/balanced function's character is greater than 99% for any of the three-qubit functions. The key ingredient 
of the read-out procedure is a direct single access to the overall set of relative phases encoded in the rovibrational 
superposition. 

IV. DISCUSSION AND CONCLUSIONS 

The present implementation of the Deutsch-Jozsa algorithm allows straightforward extension to multi-qubit binary 
functions beyond three qubits. This requires an increase in the number of rovibrational levels composing the excited 
superposition, 2" levels for n qubits. However, from the way the functions are encoded (i.e., always half the levels 
are phase and the other half 180°), an increase in n does not require an increase in the desired experimental signal- 
to-noise, i.e., there is no need to ionize more Li2 molecules. It worth noting that using laser pulses of ~1700 cm~^ 
bandwidth (i.e., ^^10-15 fs transform-limited duration) will allow the direct simultaneous coherent excitation, from 
a single selected Li2 rovibrational state of the A-state, of about 50-60 rovibrational states centered spectroscopically 
around the current excitation region of the E-state. 

The present method can be generalized for testing functions that map any even number {2m) of inputs into one 
binary (two- valued) output, and not only 2" inputs. We have chosen to work here with functions having 2" input 
values in order to create common grounds with other methods of physical implementation of quantum computation. 

In conclusion, we experimentally implemented the modified Deutsch-Jozsa algorithm for three-qubit functions using 
molecular rovibrational pure states controlled by shaped femtosecond pulses. Quantum information representing the 
function character is encoded into the relative phases of rovibrational levels that constitute the nuclear wave packet 
on the E^E+ state of the Li2 molecule. The encoded quantum information is measured directly using a weak probe 
field. This allows, in general, keeping the encoded quantum information primarily undamaged for further quantum 
operations. This is important in the context of future implementation of more complicated algorithms. 

The present task is carried out in the framework of what we consider as coherent parallel computation using quan- 
tum elements p2| . The quantum nature of the molecule is still not fully utilized for a more efficient computation. 
For example, one quantum property that can be used to advantage is the fact that the molecular eigenstates actually 
span several internal molecular degrees of freedom - electronic, vibrational and rotational; this might allow condi- 
tional manipulation and quantum projection operations. Nevertheless, this work is a very first step in exploring the 
experimental potential of using the present molecular platform for small-scale problem-specific quantum computation. 
The problem-specific approach is different from the universal quantum computation approach. It means that single 
or very few pulses will be designed specifically to induce complete unitary transformations, without the decomposi- 
tion into a circuit of basic universal gates that act on one and two qubits. For most unitary transformations, such 
decomposition is expected to produce a very inefficient quantum circuit. That is, the number of elements in the 
circuit grows exponentially with n. By small scale we currently refer to a Hilbert space composed of up to hundreds 



of quantum molecular eigenstates. This seems to be the experimentally accessible part of the very large complete 
molecular Hilbert space, which present experimental capabilities can handle. One route for larger scale computations 
is the development of much more sophisticated ultrafast molecular excitation schemes, involving, for example, well- 
controlled Rabi cycling processes. Another route to explore is a hybrid of the present technique, which is based on the 
manipulation of internal molecular degrees of freedom, with other techniques for the manipulation of external degrees 
of freedom involving, for example, trapped molecules or molecular ions |ll5|-|l7|. 

This work was supported by the Israel Science Foundation, the National Science Foundation, and the Army Research 
Ofhce. 
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FIG. 1. Schematic picture of the excitation scheme of the experiment with the relevant potential energy curves of Li2 and 
Li J [23,24]. The rovibrational structure of the E^E^ state is used for the implementation of the Deutsch-Jozsa algorithm. 
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FIG. 2. Pump-probe ionization transients originating from various tailored wave packets representing various three-qubit 
constant or balanced functions. The /i and J2 are the two constant functions, while all the other ji are balanced. The function's 
character (constant/balanced) is evaluated by measuring the signal level at a single delay time of 5 ps. 



